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A commonly used perturbative method for computing large-scale clustering of tracers of mass 
density, like galaxies, is to model the tracer density field as a Taylor series in the local smoothed 
mass density fluctuations, possibly adding a stochastic component. I suggest a set of parameter 
redefinitions, eliminating problematic perturbative correction terms, that should represent a modest 
improvement, at least, to this method. As presented here, my method can be used to compute the 
power spectrum and bispectrum to 4th order in initial density perturbations, and higher order 
extensions should be straightforward. While the model is technically unchanged at this order, 
just reparameterized, the renormalized model is more elegant, and should have better convergence 
\^ , behavior, for three reasons: First, in the usual approach the effects of beyond-linear-order bias 

■ parameters can be seen at asymptotically large scales, while after renormalization the linear model 
' is preserved in the large-scale limit, i.e., the effects of higher order bias parameters are restricted 
\ to relatively high k. Second, while the standard approach includes smoothing to suppress large 

^ . perturbative correction terms, resulting in dependence on the arbitrary cutoff scale, no cutoff- 

Q ' sensitive terms appear explicitly after my redefinitions (and, relatedly, my correction terms are less 

I sensitive to high-fc, non- linear, power). Third, the 3rd order bias parameter disappears entirely, 

, so my model has one fewer free parameter than usual (this parameter was redundant at the order 

' considered). This model predicts a small modification of the baryonic acoustic oscillation (BAO) 

^SJ ' signal, in real space, supporting the robustness of BAO as a probe of dark energy, and providing a 

, complete perturbative description over the relevant range of scales. 

(N ■ 

>■ ■ PACS numbers: 98.65.Dx, 95.35.+d, 98.80.Es, 98.80.-k 

cn : 

^ ■ 
a^ ■ 

O ' I- INTRODUCTION 

^ : 

, In the past year, significant progress has been made in using perturbation theory (PT) to calculate the large-scale 
' clustering of coUisionless mass in the Universe, to the point where essentially perfect calculations of the quasi-linear 
^-1, clustering may soon be available P, 0, Q . This motivates a fresh look at the bias models that are needed to couple 
Q ■ these calculations to most observable tracers of large-scale mass density, e.g., galaxies 0, [1], the Lya forest [1, HHIj 
' galaxy cluster/Sunyaev-Zel'dovich effect (SZ) measurements H, and possibly future 21cm surveys fl^. Baryonic 
c/3 i acoustic oscillation surveys aimed at probing dark energy [ill,[il,[ll,[ii,[ii,[ii,[l3,[il,[r9, 20, 21, 22, 23, 24, 25, 2^1, 

' in particular, fall in the range of scales where perturbation theory may be most useful. 
^ , Probably the most straightforward, commonly used, model for the bias when coupled to perturbation theory is 
to write the tracer density, 5g, as a Taylor series in the mass density fluctuations, 5 [13, H^, i.e., 5g — '^^i^iS'' 
Note that I will often refer to the tracer as galaxies (this explains the subscript g), but essentially everything in the 

■ paper could be applied to any tracer. The terms in this series generally will not decrease in size, so 5 is usually 
taken to be a smoothed version of the density field, to reduce the size of its fluctuations. When 5 is computed using 
perturbation theory for gravitational clustering [s^, HH, [H, [s^ [s^, [H, [1^ , clustering of the tracers is in principle 
fully described. [27[ showed that applying this approach to compute the power spectrum up to 4th order in the 
initial density perturbations (the lowest order that gives a correction to linear theory) leads to interesting, but not 
completely satisfactory results. They found that the 2nd order bias term (i.e., 62) produces a white (fc-independent) 
contribution to the power on large scales. Furthermore, both the 2nd and 3rd order bias parameters contribute terms 
to the effective large-scale bias, defined as the ratio of the galaxy to mass power. These terms all depend on the scale 
of the smoothing applied to the density field. Note that the most straightforward expectation is that the smoothing 
scale should be quite large, as one can see by considering the measured value of rms density fluctuations in 8h~^ Mpc 
radius spheres, erg — 0.85 [33|- This level of smoothing dramatically affects the power on relatively large scales, e.g., 
an 8h~^ Mpc radius top-hat smoothing suppresses the power at fc 0.2 /iMpc~^ by almost 50%. 

It may be intrinsically interesting, when studying galaxy formation, to determine the coefficients bi in the series 
dg = '^^biS^ and to see generally how well this model works, as a function of smoothing scale. However, I will 
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suggest a cleaner repackaging of the model, for the specific purpose of large-scale structure phenomenology, which has 
the desirable properties of preserving the form of the linear theory model at very large scales (i.e., the power at small k 
is described by a single linear theory bias parameter and a single shot-noise parameter) , and being free of any explicit 
smoothing scales. The approach I take is more or less precisely the same as the one used in classic renormalization 
of quantum field theories. When perturbation theory is pursued naively, infinities appear, but only in ways that can 
be reinterpreted as corrections to the values of parameters of the model (e.g., particle masses, analogous to our linear 
model parameters). The infinities are removed by simply redefining the free parameter in the model to be the most 
directly observable quantity (e.g., the mass you would measure by weighing the particle, or the ratio of large-scale 
g^a xy to mass power), which would be given by a series of terms if one continued to use the original parameterization 

Section |TT] of this paper gives my primary calculations, followed by some discussion in section Hill Before proceeding 
I reiterate the main goals: (1) To formulate the bias model in a way that does not require arbitrary smoothing scales, 
while simultaneously being insensitive to the small-scale, highly non-linear regime. (2) To formulate the model in 
a way that preserves the linear theory model on very large scales (low-fc), i.e., higher order bias terms only affect 
the relatively high-fc power. To be honest, I should also note here that, while I refer throughout the paper to the 
approach of [27] as the "usual" or "standard" approach, this perturbative approach to the galaxy power spectrum 
(including beyond-linear corrections) has not to my knowledge actually been used to interpret real data. However, 
between improvements in perturbation theory and the need to interpret increasingly precise observations, the time 
for this kind of approach may have arrived |26[ . 



II. CALCULATION 



We are interested in the statistics of the real-space galaxy density field, Sg (x) — pg (x) / pg — 1, with Fourier transform 
(5g(k) — j dP'x. exp(ik • x) i5g(x), where x is the comoving position. I will compute the power spectrum, defined by 
{5g{\^)5g{]^')) = (27r)35^(k + k')P<,(fc), and bispectrum, {5g{\^^)5g{\^2)5g{\^3)) = (27r)3<5^(ki + ka + k3)i?3(fci, fca, ^3). 
I am going to use perturbation theory to describe the mass density fiuctuations, 5 = 5i + 82 + 5^ + O (^f ) , where (5„ 
is of order 5" [2^, [s^, [U, IH, [H, HI, [s^ . I start by writing the galaxy density as a Taylor series in 5 [23, |2i| : 

Pa{5) =P, + P',6+ \pI 5^ + \p':,' 5' + e + {Sf) , (1) 

stopping at 3rd order because this what is needed to compute the lowest order bispectrum and the first non-linear 
correction to the power spectrum. I have added an uncorrelated noise variable e to represent shot-noise and other 
randomness in the galaxy-mass relation that appears as white noise on large scales [33|. (Generally, anything that 
affects the correlation function only at small separations will appear as a change in white noise level in the low-fc 
power spectrum. It is important to model this noise, rather than trying to avoid it by working with the correlation 
function restricted to large separations, because it contributes to the measurement errors on either statistic.) The 
variance of e is (e^) = iVg. The Taylor series coefficients and noise amplitude are effectively free parameters as long 
as we do not have a fully predictive galaxy formation model. Usually 5 is understood to be a smoothed version of 
the density field, in order to force it to be small enough to justify the Taylor expansion [13, HI]. I do not need to 
introduce this smoothing explicitly. 

As a warm-up calculation, introducing the basic idea of renormalization that I will use later, I compute the mean 
density of galaxies: 

{p,)^P^ + \pl{5'') + 0{5{) (2) 

Note that ((5'^) is not zero, but it is O {Sf) {{6} = by definition). In absence of a cutoff (smoothing scale), (5^) = 
is infinite, or at least potentially large. Equivalently, if one introduces an arbitrary cutoff, the value of cr^ will be 
sensitive to the cutoff. This is not a problem - I simply define the observed mean galaxy density to he pg = po + ^ p'^a"^ . 
I can then remove po from Eq. ([T]) in favor of pg, obtaining 

P,{^) - + P^, -5 + \(^^ {S' - a^) + \p'^' + e + O {5t) . (3) 

This calculation has not been very profound (it is so obvious that it is often done with little or no comment [13, H^), 
but it does provide a very simple example of the idea of absorbing certain kinds of bad behavior in a perturbative 
expansion into a redefinition of the parameters, in particular replacing a parameter in the original model by a more 
directly observable quantity that would otherwise be calculated as a function of the original parameters. Sufficiently 
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aggressive smoothing could of course render the 2nd term in the calculation of the mean density smaller than the 
first, but this is partly missing the point: there is no reason to allow the mean galaxy density to acquire corrections at 
each order when a simple redefinition of the original, relatively meaningless, Taylor series coefficients can guarantee 
that this fundamental observable is always directly represented by a single parameter in the model, with no need for 
smoothing. 

Note that when I say above that is "infinite or potentially large" , the uncertainty in this statement comes from the 
form of perturbation theory used to compute the mass power spectrum, e.g., standard PT [29l, Isol, [Sll. [3^ [33l. [3^. [ssl. [s^ . 
the renormalized PT of [H, or the renormalization group (hereafter RG)-improved PT of and whether and how 
much smoothing is applied. I will discuss this more below, but, in any case, it is clear that cr^ is sensitive to the details 
of one's treatment of the non-linear regime, so it is desirable to eliminate it from the final results of the calculation if 
possible. 

I now move on to fluctuations around the mean galaxy density, defining 

5^ (5) = M^Llfi ^ ci<5 + (<5^ - a^) + ic3<53 + e + O (<5t) (4) 

Pg ^ t) 

(e has been rescaled). The galaxy correlation function is 

(|Xa - Xfcl) = (5g(Xa)(5g(xf,)) = c\ {SaSb) + ^CiCg (SJ^) + ^4 ((^^(jg) - (7^) + C1C2 (SaS^) + (Caeb) + (5) 

where Sa = <5(xa). Using the fact that the 4th order terms can be treated as Gaussian at the order I am considering, 
and the linear or non-linear correlation function can be used interchangeably in these terms, 

= + ClC^a^U + \cieab + C1C2 {6a5l) + N„S^, + ... . (6) 

We see that a term has arisen, cic^a'^S^ab, that has the form of a correction to the linear theory bias, in the sense that 
there is a fc-independent factor multiplying the correlation function. It is divergent, or at least potentially large, but 
this term can be eliminated by redefining the bias to include it, similar to what I did with the mean density. I will 
not redefine the bias yet, however, because another similar term will arise from <^(5a(5^), which I have not evaluated 
yet because it requires beyond-linear perturbation theory for the density fluctuations. 
To go farther, it is simplest to move to Fourier space, where the galaxy density field is 

dg{k)^Cldk+-C2 J ^^Mk-q+gC3 J J^J^dq,6qJk-m-q2+eic + 0{5^} . (7) 

The mass density field, up to 3rd order, is 

5k = Si{k)+ ( -^-^(5i(q)(5i(k-q)J^^^(q,k-q)+ / ^^--^^-^(5i(qi)(5i(q2)5i(k-qi-q2)4^^(qi,q2,k-qi-q2) + ... 
J (27rj J (Stt) (27rj 

(8) 

where 

7(2) n 1 ^ 5 1 ki • k2 //ci fc2\ 2 /ki • kaV 

J, (kr,k2) = - + -^(^- + -j+-(^^j , (9) 

and see Eq. 11 of [27| for J^'^^ (it will not actually be used in this paper). Combining Eq. ([7]) and Eq. ([5]), we find 
the power spectrum 

Pg{k) ^ N^ + [cl+cic^a'' + —c^C2(T^]P{k) + 
^^^/(S^^(^)^(l^-^l) + 

2C1C2 j ^P(g)P(|k-q|)jf (q,k-q) + ... (10) 

Note that, if the (5's in the original Taylor series were smoothed, the term ^CiC2a^ would not take the form of a 
simple fc-independent bias, except in the limit that the smoothing scale is very small (see [40] for a different approach 
in which extreme smoothing is applied). 
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Eq. (llOp suggests the following redefinition of the linear bias parameter: 

7,2 2 I 2 , 6^ 2 

b^=c^+CiC3a + — CiC2<T . (11) 

[13] identified &f as defined here as an effective bias. The difference in approach is that I am suggesting that bi should 
now be treated as the free parameter of the model, with ci and C3 eliminated from Eq. (llOp by substitution (I write 
the result in Eq. [TBI below, after some further redefinitions to be discussed next). Note that, because the last term in 
Eq. (Uni) is aheady O (Sf) , ci in that term can be freely replaced by 61, which is equal to ci at lowest order. Note 
also that bi includes the only appearance of C3, so the third order bias has disappeared as an independent parameter. 
Eq. (jlOp contains a final divergence, or at least high-fc-sensitive term, proportional to the following: 

^P(g)P(|k-q|). (12) 

This integral diverges at high q if the asymptotic logarithmic slope of the power spectrum is rics^q — > 00) > —1.5. 
This is not the case for the standard linear theory calculation of the ACDM power spectrum; however, the fixed-point 
asymptote of the RG-corrected power spectrum of [311 is P{q) oc q^^ '^, which would give an infinite integral. As argued 
by [J, |3i 0] J the standard perturbation theory calculation based on using the usual linear theory result as the source 
for higher order terms, at all times and k, does not make much sense (and does not work very well when compared 
to numerical simulations), because at high k the initial power spectrum is quickly completely erased. Especially at 
late times, riag ~ —1.4 is much more relevant. While this is not a completely essential component of this paper, the 
beyond-linear terms in Eq. (jlOp would most naturally be evaluated using the renormalized mass power spectrum of 
Q (in fact, in the approach of Q there is really no other option). In the alternative approach of [3, Q these terms 
would probably be evaluated using the linear theory power as suppressed at high k by the renormalization of the 
propagator. 

In any case, at best the result for the term in Eq. (|12p is sensitive to the treatment of the non-linear regime, which 
we would like to avoid. Another problem with this term is that it is non-zero as A: — > 0, i.e., it is a correction to 
the power spectrum at asymptotically large scales. As discussed by [l^, this correction is fc-independent at low k, 
so it looks like a correction to the shot-noise. A basic philosophy of this paper is that we would like to preserve the 
simplest linear theory model for galaxy clustering on large scales, rather than allowing it to be corrected at each order 
in perturbation theory. Therefore, we now perform our final renormalization, absorbing this potentially divergent 
quantity, evaluated at fc = 0, into the shot-noise term, which becomes 

After this piece is subtracted, the remaining integral, written as a Taylor series around fc = 0, 

O(fc'), (14) 



(where fj, — \i - (\/kq) is clearly convergent for any reasonable power spectrum (in fact, the 0{k) term is zero after the 
fjL integration, making it even more convergent). 

The 2nd order bias is not renormalized at this order, but for notational compactness I define 



b2 = f. (15) 



The final result for the power spectrum is: 



P9 ik) = bi 



+ 4 I 7^^('7)[^(|k-q|)-P(9)] + 2 62 / ^P{q)P{\k-Ci\)jf (q,k-q) 
^ J (ztt) J (27r) 



'N . (16) 



We see that deviations from the traditional large-scale galaxy clustering model, in which the mass power is multiplied 
by a constant and white noise is added, are controlled by the single parameter, 62. The first term (linear in P{k)) in 
Eq. PB)) must be evaluated using the non-linear mass power spectrum. In standard perturbation theory, the P's in 
the extra bias terms (quadratic in P{k)) should be understood to be the linear power, while in the RG approach of 
they would be the renormalized power. In the approach of P, these P's would presumably be the initial power 
as evolved using the renormalized propagator. Note that in some technical sense there is no difference between these 
options, at the order of calculation in this paper. 
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0.1 

k [h/Mpc] 

FIG. 1: Thick lines show {Pg(k) - iV)/&?Pmass(it) for 62 = 0.6 (black), 62 = 1.3 (green), and &2 = -0.15 (red). For comparison, 
the thinner, blue, wiggly lines show the ratio of mass power in a realistic model to the power in a similar model with no baryonic 
acoustic oscillations (solid is non-linear power, dotted is linear). 

Fig. [T] shows the dependence of the power on &2- The cosmological model for this figure, plotted at 2; = 0, was 
flat with (78 = 0.85, n = 0.96, fl^n = 0.281, h = 0.71, and fib = 0.0462 with CMBfast 41] used to generate the 
transfer function. We see the desired convergence to the simple linear bias plus white noise model on large scales. 
The term oc 62 is generally negative, while the term cx 62 has the sign of 62- This means that when increasing 62 from 
zero we obtain an increase in power at first, up to 62 — 0.6, but then a decrease, especially at higher k, as the 62 
term takes over. While this bias term has a maximum, note that the simple ratio of galaxy to mass power does not, 
because extra white noise power can always be added. For increasingly negative 62 both terms are negative so the 
power decreases quickly (again, noise power can be added, so negative 62 does not automatically mean that the bias 



6 



as traditionally defined decreases). It is interesting to note that there are only hints of features related to baryonic 
acoustic osciUations (BAO) in the ratio {Pg{k) — iV)/6^Pmass(fc), at the 0.3% level for the 62 = 0.6 case. Combining 
their small size with the fact that much of the effect appears to be a modification of the amplitude of the wiggles, 
rather than their position, suggests that 1% distance measurements using BAO should not be significantly corrupted. 
Of course, these effects, the effect of broad-band /c-dependence of the bias, and the modification already present in 
the non-linear mass power spectrum, should not be ignored when fitting observations. 

In Fig. [2] I explore the sensitivity of the higher order bias terms to high-A: power. For the standard results I have 
used the RG-improved power spectrum of [3] to compute these terms, but in this figure I compare to the result when 
I use the standard linear power spectrum. We see that there is not much difference between the results for these 
cases, even though the RG power is 37% larger at fc = 0.3 /i Mpc~^, more than a factor of 3 larger at fc = 1 hMpc~^, 
and the asymptotic slopes at high k are —1.4 vs. ~ —3. This is a testament to how successful the inoculation against 
high-A; power has been (recall that the unrenormalized bias-related corrections in the RG case would be literally 
infinite, and some of the corrections using standard PT nearly so). I used 62 = 0.6 for the figure, but the results for 
larger or smaller 62 are similar. An increase in sensitivity to small scale power with increasing 62 indicates that the 
term quadratic in 62 is more sensitive than the linear term. For reference, Fig. [2] also shows the effect of Gaussian 
smoothing with rms width 2/i~^Mpc, which gives a = 1.26 for the model in the figures, i.e., not really perturbative. 
The direct power suppression for this smoothing is already substantial, reaching 30% at k — 0.3 /iMpc^^. 



A. Three-point function 



At the order in 5i considered here, the three-point function (and presumably other similar statistics that require 
only computation to this order) [H, IH, [H, [4^ 0, [13, IH, [2^ [sO, HH is described neatly by the renormalized model, 
with no non-trivial new effects. In real space it is 



((5g(xa)5g(x,,)5g(xc)) = Cf" = 4 {daSbSc) + -c?C2 {SaSb {6^ - (T^) + cyclic permutations) + 

= clCabc + ^C?C2 {^ac£.bc + Cychc) -|- ... 

Since this is already O [5f) we can freely substitute ci — > 61 to produce 



^abc 1.^ 



Cabc + (fac6c + Cychc) 



(17) 



(18) 



Deviations from the simplest model are again controlled entirely by 62. 



III. DISCUSSION 



The standard method for including biasing of tracers of mass density in perturbation theory calculations of large- 
scale clustering [27I, [1^ uses five parameters to model the power spectrum and bispectrum up to 4th order in density 
perturbations, and requires at least one arbitrary cutoff parameter (smoothing scale). In the standard approach the 
effects of the beyond-linear order bias parameters and the cutoff propagate to asymptotically large scales. (Admittedly, 
this "standard" method has not yet been used much for interpreting data, but this may change soon.) After my 
redefinitions we have only four parameters, with much more cleanly separated effects: 

standard parameters — > po, p'q, Pq, pg", A^o, and messy cutoffs 
new parameters — > pg, bi, 62, and iV . 

[27! provided the core calculations for this paper. The innovation here is the suggestion that their effective linear bias 
and shot-noise-like term should not be treated as output predictions of the model but instead should be absorbed into 
renormalized linear bias and shot-noise parameters. Significantly, these new free parameters of the model contain most 
of the high-fc (small, non-linear scale) sensitivity in the bias calculation. Rather than having a large-scale effective 
bias (ratio of tracer power to mass power) which is a function of all three bias parameters in the model, and the 
smoothing scale, I have only a single parameter for the large-scale bias. Rather than shot-noise that is a function 
of the original stochasticity, the 2nd order bias, and the smoothing scale, I have a single shot-noise parameter that 
describes all asymptotically large-scale deviation from the linear bias model. Quasi-linear deviation from the linear 
bias plus shot-noise model is described by a single parameter, with no need for smoothing. 
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k [h/Mpc] 

FIG. 2: (P9(A:)-iV)/b?Pn,ass (fe) for 62 = 0.6. The solid curve is my standard calculation, where the RG-corrected (non-linear) 
mass power spectrum has been used to compute the bias terms. To show the lack of sensitivity to small scales, for the dashed 
curve the usual linear mass power, which becomes dramatically different from the RG power with increasing k, was used to 
compute the bias terms. The dotted curve shows the effect of a 2 Mpc rms Gaussian smoothing applied to the result of the 
standard calculation. 

One might ask if it is desirable to bury time and cosmology dependence, in the form of the values of and 
J dq q^P{qf, in 61 and N. Do we not lose information this way? The answer is: not really, because Ci, C2, C3, and 
iVo are also mixed into 61 and N , and, even for fixed smoothing scale, generally none of these parameters will be 
independent of the cosmological model or of time. Furthermore, cr^ and / dq q^P{q)^ are relatively sensitive to high-fc 
power which is poorly described by perturbation theory. 

Note that simulations or other microscopic models of galaxy formation can still be used to calibrate the renor- 
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malized large-scale structure model. The way this works exemplifies the usefulness of the revised model. The most 
straightforward method for calibrating the unrenormalized model is to pick some smoothing scale and measure di- 
rectly the parameters of the original Taylor series, Eq. ([T]). If this is done very literally (e.g., by making a scatter 
plot of galaxy density vs. smoothed mass density and reading off derivatives at S — 0), one is not even guaranteed 
that the calibrated model will reproduce the mean density of galaxies in the simulations used to calibrate it (without 
the trivial rcnormalization of the first term in the Taylor series, discussed above), let alone the large-scale bias, the 
shot-noise, or really anything else about large-scale structure. (I emphasize here that I am not saying this method 
will surely fail, only that it is not guaranteed to work.) To calibrate the renormalized model one computes the mean 
density and power spectrum from the simulation, and if desired the bispectrum, and then fits for the parameters of 
the model. This guarantees that the large-scale structure in the simulation will be reproduced by the model as well as 
possible. Of course, one can be less literal and calibrate the unrenormalized model by fitting the simulated statistics, 
but this still includes the extra steps of choosing a smoothing scale and computing the simple observables (e.g., mean 
density, large scale bias) as a function of several parameters. Of course, as I discussed above, there is an additional 
problem with the literal use Eq. ^ in that any smoothing aggressive enough to render the density fluctuations small 
will also corrupt the power on interesting, perturbative scales. 

Within the perturbative approach there is a simulation-independent means of estimating the validity of the results. 
One can go through the full process of extracting the cosmological parameters of most fundamental interest from the 
data, marginalizing over the bias parameters, using the perturbative model calculated to two different orders, and 
compare the results. As presented here, the results from linear theory could be compared to the model with the first 
non-linear corrections, although computing to another order might be useful in the future (of course, higher order is 
necessary if one wants to use the bispectrum to two different non-zero orders). The level of agreement would surely 
depend on the maximum k used in the fit. This maximum k could be chosen small enough to ensure agreement, i.e., 
consistency of the perturbative model. 

There are several obvious next steps for this line of work: Extension to redshift space [H, should be straight- 
forward, following the calculations of [13]. A more general time and space-local model could be implemented, e.g., 
once one goes beyond linear theory the Taylor series for pg should probably include velocity divergence terms as well 
as density terms. As discussed by Q, it may be possible to make non-trivial predictions about the time evolution 
of bias by starting the perturbative calculation from a model for the local formation (and merging, etc.) rate of the 
tracer [Hj, [5^1 , rather than the instantaneous density of the tracer (this will probably require a rcnormalization group 
approach as in Q, rather than simple parameter redefinitions as in this paper). Extension to multiple luminosities 
or types of galaxies (or more divers tracers [H^H^, H^) can be done by writing a Taylor series like Eq. ([1]) for each. 
This will lead to extra parameters, but cross-correlation will provide extra observables to constrain them. 

Finally, I observe that, beyond just removing badly behaved perturbative terms at the present order of calculation, 
by treating the linear bias and shot-noise as free parameters we are in effect already accounting for similar terms 
that appear at any order in perturbation theory. This is probably the most compelling reason to believe that the 
approach of this paper can be expected work better than the standard approach as presented (the problem in the 
standard approach of smoothing corrupting the power on interesting scales is also significant). It may be possible to 
show that all higher order corrections to the tracer power as fc ^ can be included in these two terms, i.e., that the 
renormalized linear bias and shot-noise give a fundamentally complete descriptions of very large-scale clustering [soj . 
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